' t 

** 

* ' 

THE RESPONSE OF THE UPPER ATMOSPHERIC TEMPERATURE 
TO CHANGES IN SOLAR EUV RADIATION AND GEOMAGNETIC ACTIVITY 


by 

' S. Chandra 
B. V. Krishnamurthy* 

Laboratory for Space Sciences 
I NASA Goddard Space Flight Center 
Greenbelt, Maryland 

ABSTRACT ’ 

The upper atmospheric neutral temperature derived from 
satellite drag is analyzed with a view to determine the re- 
lative importance of solar EUV and magnetic activity on heating 
of the upper atmosphere. 

As the indices of solar EUV radiation, the radio fluxes 
obtained from the ground based observations in the frequency 
range of 200- 9400 Mc/s are used. Also for the period March 7 - 
May 15, 1962, a direct comparison is made with the EUV flux 
obtained from the OSO 1 Satellite. As the index of magnetic 
activity the planetary magnetic index K is used. 

The data have been analyzed by separating the long period 
components of 27 days and above from the short term components 
by using a 5 day running mean filter. A cross correlation study 
between the various parameters shows that the long term variations 
in temperature are strongly correlated with the similar variations 
in the radio flux in the frequency range of 1000-3750 Mc/s, 
the correlation being maximum at 1000 Mc/s. The short term 
variations in temperature, on the other hand are strongly 
correlated with those in EKp, the daily sum of Kp indices and 
not with the radio flux at any of the frequencies considered. 

The correlation between temperature and EUV flux in the short 
term components is somewhat better though not as striking as 
the long term components. 
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An empirical relation between temperature and the changes 
in radio fluxes and EK p based on the above correlation study shows 
that 

T n = 524 + 2.73 S 10>7 + 1.77 (ZK p - SK p ) 
or T n = 470 + 6.53 S 30 + 1.77 (£K p - SK p ) 

where is the nighttime minimum temperature; Is^q 7 , lS 30 

and EK are the 5 day running means of 2800 Mc/s radio flux, 

P 

1000 Mc/s radio flux and £K p respectively. 

The present investigation leads to the conclusion that 
short term variations in temperature are closely associated with 
the changes in magnetic activity represented by K and all the 
systematic variations in temperature including the so called 
semi-annual variation can be fully accounted for in terms of 
the similar variations in solar flux. 
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INTRODUCTION 

In recent years, the vast amount of observational data 
obtained from the drag analysis of a number of satellites have 
proved to be a very valuable source of information in studying 
the thermal properties of the upper atmosphere [for detailed 
references see Jacchia, 1964a; Evans, 1966; Priester et al, 1967]. 
Although the physical processes governing the temperature of 
the neutral component and in particular, the heating associated 
with geomagnetic activity, are not fully understood, the empirical 
relationships between solar decimeter flux and geomagnetic 
indices, and exospheric temperature derived from satellite drag 
observations have been remarkably successful in describing the 
thermal response of the atmosphere to changes in solar activity. 

Questions have occasionally been raised, however, about 
the suitability of a particular frequency in the radio region 
of the solar spectrum as an appropriate index of EUV radiation 
[Nicolet, 1963, Anderson, 1965]. This is because of the fact 
that the characteristics of the radio fluxes emitted from the 
sun show a marked dependence on the frequency of emission. For 
example, the amplitude of the 27 day variation is maximum at 
2800 Mc/s whereas the corresponding variation at 200 Mc/s is 
almost negligible [Smerd, 1964], The maximum variation over a 
solar cycle in the basic component is centered around 1000 Mc/s 
and falls off on either side of this frequency [Basu, 1966]. 
Unfortunately, the EUV data is available only for a limited 
period of a little over two months (March-May, 1962) and also 
for a limited range of the EUV spectrum (170A-370A) [Bourdeau 
et al , 1964] and it is not possible to make a general prediction 
for the entire solar cycle based on the data of such a short 
duration. But it is clear that even if the EUV data were avail- 
able for a much longer period, it would not be possible to 
identify a single frequency in the radio region which would fol- 
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low the EDV variation in every detail. It is, nevertheless, 
appropriate to ask, which frequency or range of frequencies 
in the solar radio spectrum can best represent the changes in 
exospheric temperature. What is the origin of the small day 
to day fluctuations in temperature which are superposed on the 
27 day and longer term variations? What causes the so called 
semi-annual variation in atmospheric temperature and density? 

To what extent can the changes in solar activity account for 
these observations? This paper is an attempt to seek answers 
to some of these questions. 

SOURCE OF OBSERVATIONAL DATA 

The temperature in the present analysis is the global 
minimum temperature (T N ) derived from the precisely reduced 
drag data of Explorer IX [Roemer, 1966] and based on Jacchia' s 
(1964b) atmospheric model. These data which approximately cover 
the period from 1961 to 1963, reflect the period of moderate 
solar activity and are particularly suited for studying the 
day-to-day variation in view of the good time resolution. In- 
cluded in the present study are also Explorer I (1958<y) data from 
1958-1963, which have been used only for the study of the long- 
term variation. The poor resolution of this data does not permit 
the study of the small day-to-day variation in a meaningful way. 

The temperature T^, has been computed from the perigee 
temperature by using the diurnal model of Jacchia (1964b) . This 
model is based on the symmetry of the diurnal bulge with respect 
to the latitude of subsolar point. There have been some discussions 
[Bourdeau et al , 1964; Jacchia and Slowey, 1966,1967] about the 
amplitude, phase and symmetry of the diurnal bulge and it should 
be admitted that there is no unique way of defining all these 
parameters with certainty. However, an uncertainty in the diurnal 
model is not likely to cause any appreciable change in the 
Explorer IX data since for the entire observation period of 
this satellite, the perigee angle with respect to the position 
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of the diurnal bulge does not exceed 90°. Thus is essentially 
independent of diurnal and spatial variation and any changes 
associated with it should be mainly due to changes in solar 
activity. 

As the indices of solar EUV radiation, the radio fluxes 
in the frequency range 200-9400 Mc/s obtained from a number 
of ground based stations and given in the following table. are used. 

TABLE 1 


Frequency of Geographic Geographic 

the radio flux Station Latitude Longitude 


200 Mc/sec 

Nera 

52 . 2°N 

5 . 1°E 

600 Mc/sec 

Humain 

50 . 2°N 

5 . 3°E 

1000 Mc/sec 

Toyokawa 

34.8°N 

137. 4°E 

2000 Mc/sec 

Toyokawa 

CO 

00 

o 

137 . 4°E 

2800 Mc/sec 

Ottawa 

45 . 4°N 

75 . 6°W 

3750 Mc/sec 

Toyokawa 

34 . 8°N 

137. 4°E 

9400 Mc/sec 

Toyokawa 

Jz; 

0 

00 

• 

CO 

137. 4°E 


Also for the period March 7 - May 15, 1962, a direct 
comparison is made with the EUV flux obtained from OSO 1 satellite 
[Bourdeau etal, 1964] . The corpuscular effects are assumed to be 
indicated by the changes in geomagnetic activities given by 
the planetary magnetic index K . 

Jr 

NEUTRAL TEMPERATURE IN RELATION TO SOLAR EUV FLUX AND MAGNETIC 
ACTIVITY 

The changes in temperature related to EUV and magnetic 
activities are studied by separating the components of 27 days 
and above from the short term fluctuations by taking five day 
running means of all the data. This method of analysis is similar 
to the one described by the authors in a recent paper (1967) . 

Plotted in Fig. 1 are the five day running means of the 
flux at 1000 Mc/s, EUV flux, T N and EK p (the daily sum of K p ) . 
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The figure shows that the 27 day variation in T N , follows closely 
the corresponding variation in EUV and 1000 Mc/s as expected where- 
as EK does not follow similar variation. What is not quite 
expected are features of Fig. 2 which shows the difference series 
of the unsmoothed and smoothed (five day running mean) data of 
the various parameters. These difference series essentially 
represent the short term variations of periods of about 5 days. 

It is seen from this figure that the fluctuations in temperature 
follow very closely with those in SK p and also to some extent 
with EUV flux, with the time lag of about 1 day. The 1000 Mc/s 
flux in general does not appear to follow similar fluctuations. 

It should be pointed out that the EUV flux used in the 
present analysis is only a fraction of the entire EUV spectrum 
responsible for the upper atmospheric heating. This may partly 
be the reason for the comparatively less striking similarity 

between T, t and EUV than between T xt and EK . The close association 
N N p 

between short term fluctuations of T„, EUV and EK raises an 

N p 

interesting question. What causes the small day-to-day fluctuations 
in the exospheric temperature? It may seem reasonable to assume 
that the small day to day fluctuations in EUV should be reflected 
in the temperature. On the other hand, it appears from Fig. 2 
that changes in the temperature associated with EK p and EUV 
are indistinguishable. If SK p is an index of solar wind velocity 
[Snyder et al, 1963] the relationship between EK and EUV if 

r 

any, should reveal a time lag of about 2 days or more. Such a 
time difference is not always apparent from Fig. 2. Could it be 
then that some of the small day-to-day fluctuations in EK are 

r 

partly the result of the atmospheric circulation in the lower 
atmosphere due to changes in temperature which in turn are caused by 
changes in EUV flux. This possibility was suggested by Newell (1966). 
It would be difficult to arrive at any definite conclusion about 
this hypothesis from the present analysis in view of the limited 
amount of data and the difficulty of separating the contributions 
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from terrestrial and extraterrestrial causes. It nevertheless is rea- 
sonable to conclude that EUV plays an important role in affecting 
the upper atmospheric temperature both for short and long term 
variations and not withstanding the lack of full interpretation 
°f Kp, it is a good index of short term variations in temperature. 

CORRELATION BETWEEN T M , £K AND SOLAR FLUXES AT DIFFERENT 
FREQUENCIES w P 

In spite of the close association between temperature and 
EUV as revealed from the preceding analysis, an empirical relation 
cannot be established between the two parameters in view of 
the limited amount of data available for EUV. Instead, it is 
possible to investigate the relationship between temperature 
and solar radio fluxes at different frequencies since continuous 
data are available for both radio fluxes and temperature. 

Using the Explorer IX neutral temperature data from the 
period Nov. 1961 to May 1963, correlations between temperature 
and radio fluxes at different frequencies are obtained giving 
a time difference of -5 to +5 days between the parameters in 
steps of one day after smoothing the data by taking five day 
running means. The maximum correlation which shows a time lag 
of one day between T^ and the fluxes is plotted as a function 
of frequency in Fig. 3. It is seen from this figure that the 
correlation rapidly increases with frequency and attains a 
maximum value at about 1000 Mc/sec. Afterwards it gradually 
falls off at higher frequencies. Since the correlation at 
2800 Mc/sec is only slightly less than at 1000 Mc/sec, it seems 
both these fluxes can be used for obtaining an empirical relation 
between temperature and the fluxes. In fact, any frequency in 
the range 1000-3750 Mc/s may be used without an appreciable error. 

The correlation between the short term fluctuations (differences 
between unsmoothed and 5 day running mean series) in T^ and 
any of the radio fluxes is negligibly small. On the other hand, 
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the correlation between the short term fluctuations in T, t and 
EK p is very high as was shown in an earlier paper [Chandra and 
Krishnamurthy , 1967]. Based on the short and long term 
correlations between the various parameters, the following 
empirical relations are obtained (using the linear regression 
formula) . 

V = 524 + 2.73 S 1() ? + 1.77 (SK p -EK p ) (1) 

or t n ’ - 470 +6.53 S 30 +1.77 (ZK p -SK p ) (2) 

where "S 30 , 7 , and ZK p are 5 day running means of 2800 and 

1000 Mc/sec flux and £K respectively. In view of the one day 

hr 

time difference between and the various parameters the T^' 
calculated by using the above relations, should be assigned to 
the next day. We wish to emphasize here that EK p is being used 
only as an index to represent the short term variations in 
whatever be its real physical significance. 

COMPARISON BETWEEN TEMPERATURES OBTAINED FROM THE EMPIRICAL 
RELATION AND FROM SATELLITE DRAG DATA 

The temperatures calculated from the empirical relations 
given in the preceding section are now compared with the 
Explorer IX temperature data for the period November 1961 to 
May 1963. Although the same Explorer IX data has been used 
to obtain the empirical relations, the comparison would still 
be meaningful as the correlation coefficients from which the 
empirical relations are obtained, are less than unity. 

Plotted in Fig. 4 are the differences T^-T^' where Tjj 
are the Explorer IX temperature and T^’ are those obtained from 
equations (1 and 2). It can be seen from this figure that the 
residual Tj^-T^' , in general, varies between -50° to +50°K and does 
not . show any systematic variation. This places a fair amount 
of confidence in equs. ( 1 and 2) in their ability to represent 
all the systematic variations. 
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As an additional test of equs. (1) and (2), comparison 
is also made between the temperature computed from these equations and 
Explorer IX temperature data for the period Feb. 1961 to Oct. 1961 
which were not included in formulating the empirical relations. 

This period is particularly interesting for comparison because 
of the large changes in temperature associated with geomagnetic 
activities. It was found that the comparison is extremely good 
even to the finer details. This is illustrated in Fig. 5 which 
shows the plot of T^' computed from 1000 Mc/s flux and the 
Explorer IX temperature data as given by Jacchia and Slowey (1965) . 
Because of the difference of model, a correction of 5% is applied 
to this data to make it compatible with the Jacchia' s atmospheric 
model (1964b) which was used in obtaining Explorer IX temperature 
for the period Nov. 1961-May 1963. An examination of Fig. 5 
clearly shows the remarkable similarity between the two temperatures. 
Even during the periods of high magnetic activity, the fluctuations 
in both Tjj and T^' , in general, follow each other very closely. 

Thus the variations in temperature associated with magnetic 
activity are well represented by the K index both during the 

r 

quiet and disturbed periods. 

From the comparisons made in this section it is evident 
that the empirical relations presented in the previous section 
are quite adequate to represent the global minimum temperature 
of the neutral atmosphere above thermopause. 

SEMI-ANNUAL VARIATION IN TEMPERATURE 

It was shown in the previous section that the empirical 
relation based on solar flux and K represents all the short 

IT 

and long term variations in temperature and the residual T^-T^' 
does not show any systematic variation. It has been reported 
by a number of workers [Paelzold and Zschorner, 1961, Jacchia 
1965, King-Hele, 1966, Cook; and Scott , 1966;; Cook,r.T967] 
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that the atmospheric density and temperature after correcting 
to a standard solar activity show a semi-annual variation. The 
physical explanation for this effect has hitherto remained a 
mystery though a number of mechanism have been proposed [Paetzold 
and Zschorner 1961, Priester and Cattani 1962, Jacchia 1965, 

Anderson 1966, Cook 1967], Since the residual T N -T^' in Fig. 4. 
does not show any systematic semi-annual variation, it appears 
that the so called semi-annual variation could be accounted 
for in terms of the solar flux variations. This point is in- 
vestigated in further detail in the following. 

Using again the Explorer IX T N data, the solar fluxes at 2800Mc/s 

and 1000 Mc/s and £K for the period November 1961 to May 1963, 

P 

the 27-day and shorter components are smoothed out by taking 
27-day running means. The resulting smoothed data essentially 
contain components of periodicities above 4 months. Plotted 
in Fig. 6 are the 27 day running means of all the four parameters. 

It can be seen from this figure that there is a remarkable 
similarity between the variations of Tjj and those of S^q and 
S^ 0 rj. All these three parameters show broad peaks around the 
months of March and October. The variations in EK do not 

r 

follow those in T, t , or S OA or S, n _. Quantitative estimates 
of the degree of association between the long term variations 
in the various parameters under consideration are obtained by 
evaluating the correlation coefficients, shown in Table 2. 


Parameters 


T N vs S 30 
T N vs S 10.7 
T N ys SK p 


Correlation 

Coefficient 

0.79 

0.76 

0.16 


TABLE 2 
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It can be seen from Table 2 that shows a high degree 
of correlation with both S^q and S^q ^ whereas the correlation 
with SKp is poor. To investigate these long term variations 
over a longer period data from lSSSa satellite which are avail- 
able from 1958 to 1963, are analyzed along with the 2800 Mc/sec 
solar flux data from the same period. Components of 27 days and 
shorter periodicities are eliminated again by taking 27 days 
running means of the data. The resulting smoothed data of T^ 
and ^ are plotted for each year separately as shown in 

Fig. 7. Also shown in the same figure are the smoothed data of 

the Zurich sunspot number (R ) for each year separately. 

z 

Fig. 7 reveals a remarkable similarity between the variations 
in T n , S^ 0 ^ and R z . It is interesting to note that the semi- 
annual variation in T^, ^ , and R z is not a persistent feature 

throughout the period from 1958 to 1963. Marked semi-annual 
variation is shown only during the years 1958 and 1962 and 
indicated in 1963. During the period 1959-1961, all the three 
parameters show oscillatory behavior and there is no clear trend 
of semi-annual variation though occasionally peaks around April 
and October could be detected among other peaks. In making 
the comparison between temperature and 10.7 cm flux, it should 
be remembered that the former is derived from the perigee temperature 
using an empirical diurnal model. Some of the phase differences 
between the corresponding peaks may be attributed to the uncertainties 
in the diurnal model . 

The similarity between the long term variations in T^, S^q ^ , 
and R strengthens the view that they are not spurious. It is 
evident from the present analysis that the long-term variations 
in temperature can be fully accounted for by the similar variations 
present in the solar flux and it is not necessary to invoke any 
other mechanism to explain the long-term variations in temperature. 

The physical mechanism giving rise to these fluctuations in solar 
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fluxes is not within the scope of this paper and will be discussed 
elsewhere. From the aeronomical point of view there appears 
to be no difficulty in understanding the so called semi-annual 
variation . 

SUMMARY AND CONCLUSIONS 

The present investigation of the variations in neutral 
temperature above thermopause in relation to solar and geo- 
magnetic activities has led to the following conclusions. 

1. The short-term fluctuations in temperature are associated 

with similar fluctuations in both EUV and EK . The correlation 

P 

between temperature and EK is however greater than that between 

P 

temperature and EUV. 

2. The variation in temperature on a long-term (27 days 
and above) basis shows a good correlation with the similar 
variations in solar flux in the frequency range of 1000-3750 Mc/s , 
with maximum correlation at 1000 Mc/s. 

3. The following empirical relations seem to account for 
all the systematic variations in temperature both during quiet and 
disturbed conditions. 

T n = 524 + 2.73 S 1Q ? + 1.77 (SK p -EK p ) 

or T n = 470 + 6.53 S 30 +1.77 (SK p -EK p ) 

4. The so-called semi annual variation in the upper 
atmospheric temperature and density is not a persistent feature 
throughout the solar cycle and may be accounted for in terms 

of similar variations in solar activities without the necessity 
of invoking any other involved mechanism. 
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FIGURE CAPTIONS 

1. Daily plots of 5 day running means of 2K p , T N , EUV and S 3Q . 

2. Daily plots of the differences between original and 5 day 
running mean data of EK p , T N , EUV and S 

3. Correlation of temperature with radio fluxes at different 
frequencies . 

4. Plots of the residual T xt -T ' . 

JN N 

5. Comparison of T^' computed from 30cm flux model with Explorer IX 
temperature for the period February 1961 to October 1961. 

6. Daily plot of 27 day running means of S^q, T^> ^ and EK p . 

7. Daily plot of 27 day running means of Tjj, ^ and R z . 
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